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■ Abstract 

In this paper we construct a family of exact strong solutions to the two- 
dimensional incompressible liquid crystal equations with finite energy. The ini- 
tial velocity is chosen to be rotationally symmetric and the image of the initial 



in 
o 

(N 



orientation of the liquid crystal is a non-trivial curve on the unit sphere. It turns 
out that this family of initial data evolves globally in time by liquid crystal flow 
and may shrink to a single point as time goes to infinity. 



<N . 1 Introduction 

> 

We consider the following hydrodynamic system modelling the flow of liquid crystal 
materials in two dimensions (see, for instance, [3], El [7]): 



u t + u ■ Vu + Vp = Au - V ■ (Vd <g> Vd), 
{ d t + u- Vd = Ad+ |Vd| 2 d, (1.1) 
^. V-u = 0, |d| = l, 

c5 ■ where u is the velocity field, p is the scalar pressure and d is the unit-vector on the 
sphere S 2 C M 3 representing the macroscopic molecular orientation of the liquid crystal 
materials. Here the ith component of V ■ (Vd £g> Vd) is given by Vj(Vid ■ Vjd). For 
simplicity, we have set all the positive constants in the system to be one. We are 
interested in the Cauchy problem of (II. ip with the initial data 

u(0, x) — Uo(a;), d(0, x) = d (x). 
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The above system (II. ip is a simplified version of the Ericksen-Leslie model for 
the hydrodynamics of nematic liquid crystals [3j |6]. The mathematical analysis of 
the liquid crystal flow was initiated by Lin and Liu in [TJ E] . The existence of weak 
solutions in two dimensions is obtained in [10] where the authors also showed that 
there are at most finitely many time singularities for their weak solutions (see also jl]). 
The uniqueness of weak solutions in two dimensions was studied in [Hi [13] . In [12] , 
the global existence of strong solutions is proved for initial data with sufficiently small 
norms of u and Vd in BMO -1 . See also [9] for a small data global existence result 
in 3D. The global existence of large solutions is obtained in two dimensions in [2j [10] 
under the assumption that the image of do is contained in a half sphere (see also a 
recent new proof in [5] ) . 

In this paper, we are concerned with a family of exact strong solutions with large 
initial data, which are global in time. We succeed in constructing them by choosing the 
rotationally symmetric initial velocity fields and symmetric initial orientations. Thus 
the image of the orientation must be a curve on § 2 which could be non-trivial. We also 
show that the image curve of d on § 2 may finally shrink to a single point as time goes 
to infinity. 

The main result of the paper is the following Theorem: 

Theorem 1.1. Let Si G (0,7r/2) and Uo,ipo be two given functions of r G K + which 
satisfy 

/*0O /"OO 

/ ul(r)r dr < oo, / (i/j' (r)) 2 r dr < oo Si < ^o( r ) < tt — 5\. (1.2) 
Jo Jo 

Then there exists a smooth function $ : (Si, ti — Si) — >■ R such that, for the rotationally 
symmetric initial velocity field 

u (x) = u (\x\)e e , e e = (-^:,^T) T (1.3) 

|x| \x\ 

and the symmetric initial orientation unit vector 

/sin^o(M)cos M\ x \)\ 
d (x) = smifj (\x\)sm(f) (\x\) , O = $(^o), (1.4) 

\ COS^o(M) / 

the solution to the liquid crystal equations ( II. ip with initial data ( ll.3p - (11.4p is global 
and can be exactly solved. 

Remark 1.1. The solution we obtained in Theorem 11.1 1 becomes smooth instantaneously 
for t > 0. By the weak-strong uniqueness (see, for instance, [11] and [13]), such solution 
is also unique. 

We point out here that the L 2 norms of u and Vd may be large in our theorem. 
Hence the global well-posedness of such solutions is not a consequence of [12]. Since 
the image of do is not contained in any half sphere on S 2 (cf. ( 13. 6p ). the existence of 
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our solutions cannot be deduced from [21 QUI E] either. The main ingredient of this 
paper is that we can give an explicit expression of the solution to the liquid crystal 
equations ( II. ip with initial data fll.3p -( TT~4]) . For the proof, we first parameterize the 
curve, and solve a nonlinear ordinary differential equation (12. 12}) to reduce the problem 
to a nonlinear heat equation. We then construct a family of exact solutions by solving 
the nonlinear heat equation using a Hopf-Cole type transformation. We also show 
that the image curve of d has non-increasing arc length along the liquid crystal flow. 
However, the arc length does not always decreasing in general. We will give an example 
of global solutions which does not shrink in time at the end of this paper. 

Before ending this introduction, let us mention a related result on harmonic map 
heat flow in [TJ where finite-time singularities are shown for a class of initial data in 
the form 

(X\r~ x sin ipo(r)\ 
x 2 r~ l sin?/>o(V) 
cos^o(^) / 

with 

^o(O) = 0, i> (R) > vr for some R > 0. 

This also serves as an example of finite-time singularities for the liquid crystal system 
( II. ip in the special case when u = 0. Whether there are finite-time singularities for the 
incompressible liquid crystal flow in two dimensions with finite energy and non-trivial 
velocity remains a very interesting open question. 

The paper is organized as follows: In Section 2 we will parameterize the orientation 
unit vector d and rewrite the liquid crystal equations ( II. ip in terms of <ft and r/>. Then 
in Section 3 we solve ( II .ip by giving an exactly expression of u and d, and finish the 
proof of Theorem 11.11 

2 Parameterization of the Liquid Crystal Equations 

In this section we parameterize the liquid crystal equations. Let e r = x/\x\. It is clear 
that the gradient operator can be expressed as 

V = e T d r + r~ x e e d dl d e = x±d 2 - x 2 di. 
We look for solution of (u, d,p), where u is rotational and p, d are symmetric, i.e., 

d = d(£, r), u = u(t,r)eg, p = p(t,r). 
We can rewrite (11.11) into the following system 

u t = u rr + -u r — -^u 
< p r = i M 2 -|(|d r | 2 ) r -Ad-d r (2.1) 
d t = d rr + ±d r + |d r | 2 d 
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Since d £ § 2 , we may assume 



''sin ip cos N 

sin?/' sin ] . (2.2) 

COS?/' 



A simple computation gives 





( cos ip COS 0^ 




f — sin-0 sin<A 




cos sin (f) 




sin -0 cos 




{ — sin ip J 




V o / 



(2.3) 



and 

(— cos^sukA /— sin?/' cos( 

cos^cos</> , d H = -sinT/'sin^ | . (2.4) 
) V o 

Clearly, from (Q, ((23]), and flU}, we have 

d^, • d = 0, • d^, = 1, d i) d (f) = 0, 
d^ • dw, = 0, d^ • d^ = 0, d^- d^ = - cos ip sin ■0, 

and 



d^ • d = 0, d^ • d^ = 0, d^ • d^ = sin 2 ip, 
d^-d^ = 0, d^ • d^ = sin^cosV', d^ • d^ = 0. 



By the chain rule, 



(2.5) 



(2.6) 



d t = d^ip t + d<p(j) t , 
d, = d^ipr + d o,.. 

(Vv) 2 + 2d 1 p ( f ) 1p r (f) r + d^((/) r ) 2 + d^1p rr + d^rr- 

Combining the above equalities with the third equation of ( 12. ip . we obtain 
d^ip t + d^ t = d^(ip r ) 2 + 2d^ip r (j) r + d^(0 r ) 2 + d^ip rr 

+d^<p rr + -(d^ip r + d^ r ) + d{if$ + sin 2 ^ 2 .)- (2.7) 

We dot (12. 7p with d^ and d^, respectively and use (12. 5 p and (12. 6 p to deduce the following 
two equations 

ip t = ip rr H — ip r — cos ip sin ipcpl, (2.8) 
r 

6 t sin 2 = (<p rr H — r ) sin 2 ip + 2 cos ^ sin ipip r (p r . (2.9) 
r 

Note that in (12.7p the terms in the d-direction cancel each other. 
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We assume that the image of d is a curve on § 2 . In other words, <p and ip are not 
independent. Take <p = $(V0- Then it follows from f |2.8j) and (12. 9p that 

ip t = i/j rr + -%p r - ($') 2 cos ip sin (2.10) 
r 

$Vt sin 2 ^ = $'(VVr + ~Vv) sin 2 ip + &'(ip r ) 2 sin 2 ip + 2$' cos ip sin ip(ip r ) 2 ■ 

r 

By substituting the expression of ip t in (I2.10p into the last equation, we then deduce 

(&f cos ip sin 3 ip(ip r ) 2 + $"sin 2 ^(Vv) 2 + 2$' cos^ sin^(^ r ) 2 = 0. (2.11) 
A sufficient condition for (12. lip is the following ODE 

$" sin 2 ip + 2$' cos ip sin ip + ($') 3 cos ip sin 3 ip = 0. (2.12) 

3 Proof of Theorem 11.11 

In this section, we complete the proof of Theorem II .11 by explicitly solving the nonlinear 
ODE (I2.12p and a nonlinear heat equation. Note that (12.121) is equivalent to 

(sin 2 ip&)' + cos ip sin 3 ip(&) 3 = 0. 

Thus, for ip 7^ 0, tt, we have 

(sin 2 ip<&')' cosip 
(sin 2 ip&') 3 sin 3 ?/'' 

and then 

1 + i> 



(sin 2 ip&') 2 sin 2 -?/; 
for some constant (3 > 1. Consequently, 

(^) 2 = (R ■ S , 1 ■ 5 7 - (3-1) 

(p sin -0 — 1) sin 

Here we require 

/3sin 2 ^>l. (3.2) 
Recall the condition on the initial data 

tpo £ (Si, ir — Si). (3.3) 

We fix a /3 > 1/ sin 2 5 X . Plugging the expression (13.11) into (12.101) . we get 

1 COS Ip 2 
A = Vrr + ~1pr ~ /l0 • 2 / TT^ F^r) ( 3 - 4 ) 
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which is a nonlinear heat equation in 2D. This motivates us to use the idea of the 
Hopf-Cole transformation to reduce f l3.4[) to a linear equation. Let F = F(ip) be a 
function to be chosen. Clearly, 

F t = Ftyt, F rr + ^F r = F'(^p rr + -V) + F> 2 . 
Then F satisfies the linear heat equation 

F t = F rr + ~F r (3.5) 

r 

provided that 

F" + F'- = 0. 

(/Ssin if) — 1) sin-0 

We solve the above ODE to get 

sin^ 



(/^sin^-l) 1 / 2 
and 



*W = / ^ (/3si j"%i/ 2 # = ^72 ( Cl arccos (/^T cos ^) + 672 

Here Ci and C 2 are arbitrary constants. For our purpose, it is convenient to take 
d = (3 1 / 2 and C 2 = 0. Then 



F(ip) = arccos 



Note that F is well defined for any ip G (5i,n — Si), and it is a strictly increasing 
continuous function in the same interval. By ( 13. 31) . we have 



FOq) e (5 2 ,vr-5 2 ), 5 2 = arccos (\J cos fr.) £ (0,vr/2). 
From (13.51). we have 



F^(t,-)) = v(t,-)*F^ ), 

where T(t, •) is the 2D heat kernel. Hence, 

F(lf>(t, •)) £ (fc, 7T - (Ja), •) e 7T - <J X ), 

which implies that the condition (13. 2p is satisfied for any t > if it is satisfied at t — 0. 
Therefore, the solution ip exists globally. Furthermore, by the maximum principle, the 
length of the image curve of d is non- increasing. 
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It follows from (13. 1 j) that 



(/3sin ^ — 1 ) - 1 / 2 sin t/> 



Here we took the positive branch of the square root. Therefore, by choosing a suitable 
constant, 

$ ^ = / To - ~~2 lU/2 • ^- 

A/2 (P sin s - lj-V^sms 



Observe that 



lim $(*0) -> -oo, lim $(^) ->■ +oo as <5i ->■ 0. (3.6) 

if)— >Si Ip—tTT—Sl 



We now treat the first two equations in ( 12.11) . In the equation of u, we divide both 
sides by r and get 

/,M +-(-). (3.7) 



. r / t \r / rr r \r / r 



We note that d rr + is the Laplace operator in R for radially symmetric functions. 
Moreover, as a radial function in R 4 , Uq/t G L 2 (R 4 ) since, by (II. 2p . 

POO 

u l{\v\)\v\~ 2 dv = c ul(r)rdr < oo. 



Therefore, due to (13. 7j) we get 

u(t,r) = r e « u (t, | S/ 1 ) 1 2/ 1 1 dy, x = r(l, 0, 0, 0), 

J R 4 (4vrt)^ 

and, as a radial function in R 4 , u(t, -)/r G L 2 (R 4 ) for any £ > 0. This further implies 
that u(t, •) G L 2 (R 2 ). Since u = (—x 2 ,Xi)u/r, we see that u is a smooth function. 
Finally, from the second equation of (12. ip . we get 

p r = -u 2 — d r ■ d rr — (d rr H — d r ) • d r = —u 2 — 2d r ■ d rr . d r ■ d, r . 

ry ry ry ry 

Therefore, 

p(t,r) = j (^u 2 (t,s)-^\d r (t,s)\ 2 ^j ds-\d r (t,r)\ 2 . 

Note that by the radial symmetry condition, the integral above is convergence. This 
completes the proof of Theorem 11.11 

Remark 3.1. In the special case that the limit tpo exists as r 4 oo, from the proof 
above it is easily seen that i[>(t,-) converges to this limit uniformly in R 2 as t goes to 
infinity. Indeed, F(ijj ) is a bounded function and has a limit F as r — > oo. By the 
simple property of solutions to the heat equation, we have 

lim F(ip(t, ■)) -»■ F 

t— ¥OD 



uniformly in M 2 , which implies that 

lim ip(t, •)—>■■?/>:= arccos 



t— >co 



'0 



(3 



cos F 



uniformly in IR 2 . In this case, the image curve shrinks to a point on § 2 uniformly along 
the liquid crystal flow. 

On the other hand, one can find certain initial data such that the image curve on 
§ 2 does not shrink at all. This can be seen from the one-to-one correspondence of F 
and ip and the following example. Let vq be a radial function in R 2 defined as 



i 



v (r) = on [0,e], v' (r) = < 



i _ 

r (6fc+2) 3 -(6fc+l) 3 
\3 



on e 



(6fc+l) 3 (6fc+2) 3 l 



^ r (6fc+5) 3 -(6fc+4) 3 

It is easy to check that Vv G L 2 (R 2 ), v e [0, 1], and 

1 On [ e (6^2) 3 )e (6fe+4) 3 - 



On (e( 6fc + 2 ) 3 , e( 6fc+4 ) ] U ( e (6fe+5) 3 ; e (6/c+7) 3 ] ^ ken. 
On ( e (6fc+4)3 ? e (6fc+5) 3 ] 



Vo(r) = 



On [ e (6fe+5) 3 )e (6fe+7) 3 



,JfcGN. 



Now take tk = e 2( - 6k+3 ^ , tk = e 2 ^ k+6 ^ , and let w be the solution to the heat equation 
with initial data v . Since 



g 2(6fc+2) 3 e 2(6fc+4) 3 

lim : = 0, lim 



fc— s-oo ti 



k— >oo 1 1 



OO, 



we get 



lim 

fe— >co 



r(t fc , x) (ix = i. 

| x | e ( e (6fe+2) 3 ie (6fe+4) 3 ) 



Therefore, 
Similarly, we have 
Thus, the range of v does not shrink in time. 



lim v (tk, 0) = 1. 

k^oo 



lim v(t k ,0) = 0. 

k^oo 
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